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The quantitative analysis of the "catastrophic relaxation" of the coherent spin precession in 3 He-B is 
presented. This phenomenon has been observed below the temperature about 0.5 T c as an abrupt shortening 
of the induction signal decay. It is explained in terms of the decay instability of homogeneous transverse 
NMR mode into spin waves of the longitudinal NMR. Recently the cross interaction amplitude between the 
two modes has been calculated by Sourovtsev and Fomin for the so-called Brinkman-Smith configuration, 
i.e. for the orientation of the orbital momentum of Cooper pairs along the magnetic field, L || H. In their 
treatment, the interaction is caused by the anisotropy of the speed of the spin waves. We found that in the 
more general case of the non-parallel orientation of L corresponding to the typical conditions of experiment, 
the spin-orbital interaction provides the additional interaction between the modes. By analyzing experimental 
data we are able to distinguish which contribution is dominating in different regimes. 



PACS: 67.57.Lm, 76.50.+g 

1. INTRODUCTION 

For magnetically ordered systems the instability of 
homogeneous precession is a well known phenomenon. 
Suhl [2] explained it in terms of parametric instability 
of the mode of precession with respect to excitations of 
pairs of spin waves satisfying the condition of resonance: 



nuj L = w s (k) -I- w s (-k) 



(1) 



where wl is the precession frequency and n is integer 
(see also book j^] ). Parametric instability (see, e.g. 
textbook by Landay and Lifshts, 0]) is a particular 
case of the decay instability, which is very general and 
well known in physics of nonlinear wave phenomenon, 
that includes induced light scattering (photon decaying 
into photon and phonon or spin wave) decay of Leng- 
muir waves into Lengmuir and ion-sound waves in non- 
isotherm plasmas, decay of capillary wave into capillary 
and gravity wave on the water surface, etc. In quan- 
tum solids and liquids the decay instability has been 
observed in anti-ferromagnetic solid 3 He 5 , and has 
been predicted for superfluid liquid 3 He-A [S] where it 
has been observed later 171. 
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Unique feature of superfluid 3 He-B is the phe- 
nomenon self sustained and long-lived precession with 
the coherent phase across the whole precessing do- 
main, the so-called Homogeneously Precessing Domain 
(HPD) 8 . What is the really amazing is a huge preces- 
sion angle up to 104° in the presence of lower- frequency 
spin waves (longitudinal NMR), for which the resonance 
conditions arc definitely satisfied. In general nonlin- 
ear media the life time of this kind of excitation (with 
dimensionless amplitude of order unity) should be about 
their period. The main question is why the HPD is so 
long-living excitation? The explanation is unique sym- 
metry of the leading Zeeman interaction which therefore 
does not contribute to interaction amplitude of the de- 
cay processes Q [that we denote as V(k)\. Moreover, 
even subleading spin-orbit interaction in the Brinkman- 
Smith configuration (with L || H) does not contribute 
to the interaction amplitude V(k) due to the symmetry 
constraint. 

Nevertheless the abrupt instability of HPD and of 
the homogeneous precession in general, called the catas- 
trophic relaxation, has been observed in 3 He-B below 
~ (0.4 — 0.5)T C [§]. The main physical question here is 
to clarify what is the origin of contribution to interaction 
amplitude V(k) that is responsible for the catastrophic 
relaxation in general and in the experimental conditions 
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P| in particular. First reasonable step in this direction 
was recently made in pQ who considered simple L || H 
configuration and found non-vanishing contribution to 
V(k), originated from the dependence of the spin- wave 
velocity on the direction of propagation. As we will 
show below this contribution hardly can be considered 
as the main one in typical experimental conditions 

In our Letter we found another contribution to the 
decay-interaction amplitude V(k) [denoted below as 
V BISV (fc)] that allows one to rationalize main features of 
observed catastrophic relaxation. Our point is that un- 
der conditions of the experiment [H], the boundary con- 
ditions on the wall of container induce the texture of 
the order parameter in which the orbital vector L devi- 
ates from its symmetric orientation along the magnetic 
field H in the most of the container volume. The sym- 
metry of the spin-orbit interaction is violated providing 
the additional term in the interaction V(k) between the 
modes, which is dominating in typical experiments with 
the catastrophic relaxation pi llUlITT] . 



2. GENERAL PROCESSING STATES 

The homogeneous precession of magnetization in 
3 He-B has been analyzed by Brinkman and Smith 
and in a great detail by Fomin f° r the ideal case 
of the symmetric orientation of the orbital momentum 
L || H. To discuss the real experiments in which 
the orbital vector L is deflected due to the boundary 
conditions and forms the texture, the computer sim- 
ulations on the basis of the full set of Leggett-Takagi 
equations for spin dynamics have been employed, us- 
ing the program elaborated by Golo ^3] f° r the one 
dimensional texture. The instability of homogeneous 
precession was demonstrated: it was found that at the 
temperatures, corresponding to the experimentally ob- 
served catastrophic relaxation, some mode with fixed 
non-zero wave vector k starts to grow exponentially. It 
was also found JS] that the magnetic field dependence 
of the onset of instability is in quantitative agreement 
with the field dependence of catastrophic relaxation ob- 
served in ^B]. This demonstrates that the catastrophic 
relaxation is indeed in the frame of the Leggett-Takagi 
equations, but for the real understanding of this phe- 
nomenon we must identify the principle mechanism of 
this calculated instability. For that we modified the the- 
ory for the general case of arbitrary orientation of the 
orbital vector L. 



21 Symmetry of precessing states 

Let us consider the general homogeneous free pre- 
cession in external magnetic field H. In liquid 3 He the 
spin-orbit (dipole-dipole) interaction is weak. If it is 
neglected, we can apply the powerful Larmor theorem, 
according to which, in the spin-space coordinate frame 
rotating with the Larmor frequency the effect of mag- 
netic field on spins of the 3 He atoms is completely com- 
pensated. Since the magnetic field becomes irrelevant, 
the symmetry group of the physical laws in the precess- 
ing frame is 

G = SO% x SOl , (2) 

where SO^ is the group of orbital rotations in the 
laboratory frame; and SO§ is the group of spin ro- 
tations in the rotating frame whose elements g(t) are 
constructed from the elements g of conventional spin 
rotations in the laboratory frame: 

g(t) = 0- 1 (i,iu L t) gO(z,w L i) . (3) 

Here the matrix O a p(z,u)t) describes the transforma- 
tion from the laboratory frame into the rotating frame 
- this is the rotation about the magnetic field axis z by 
angle uj^t. Now we can find all the degenerate coherent 
states of the Larmor precession applying the symmetry 
group G to the simplest equilibrium state of the given 
superftuid phase: A = O^R^OAW) (rC 2 ))- 1 , where 
R*- 1 ) is the arbitrary matrix describing spin rotations 
in the precessing frame and R*- 2 -* is another arbitrary 
matrix which describes the orbital rotations in the lab- 
oratory frame. In case of 3 He-B this state corresponds 
to the state of Cooper pairs with L = S = 1 and the 
total angular momentum J = |17|: 

A<$=A B 6 ai . (4) 

The action of elements of the group G on this stationary 
state leads to the following general precession of 3 He-B 
with the Larmor frequency (if the spin-orbit interaction 
is neglected): 

A ai (t) = AR ai (t) , (5) 

R ai (t) = O a ^z,-ut)R^O^(z,cot)(R^)-l . (6) 

The matrix R^ 1 ' determines the direction of spin density 
in the precessing frame: 

S a = X R§Hf, , (7) 

where \ is the spin susceptibility of 3 He-B. This corre- 
sponds to the precession of spin with the tipping angle 
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cos (3\ = Rzz ■ The matrix R/ 2 ) determines the direction 
of orbital momentum density in the laboratory frame: 



Li — —R ai (t)S a (t) 

(2) 

with the tipping angle cos/3 2 = R Z z ■ 

22 Spin-orbit interaction as perturbation 

The spin-orbit interaction in 3 He-B is 



(8) 
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cos6»(t) + -j , 



(9) 

where £l L is the so called Leggett frequency - the fre- 
quency of the longitudinal NMR; 6 is the angle of rota- 
tion in the parametrization of the matrix R a i in terms 
of the angle and axis of rotation ; we use the system 
of units in which the gyromagnetic ratio 7 for the 3 He 
atom is 1, hence the magnetic field and the frequency 
will have same physical dimension. 

In the general state of the Larmor precession (JjJJ, 
the spin-orbit interaction contains the time independent 
part and rapidly oscillating terms with frequencies u>l, 
2u>l, 3wl and 4ujl' 



(10) 



The time-independent part - the average over fast os- 
cillations - is 

Fo = ^X&sl -\ + \& + + 0cos 7 ) 2 + 

i(l- S ) 2 (l-0 2 + (l-^ 2 )(l-^)(l + cos 7 )] • (11) 

Here s = cos fa and I — cos fa are z projections of unit 
vectors s = S/S and 1 = — L/L; and 7 is another free 
parameter of the general precession. Altogether the free 
precession is characterized by 5 independent parameters 
coming from two matrices and R/ 2 ' ^Hj : two an- 
gles of spin S, two angles of the orbital momentum L, 
and the relative rotation of matrices by angle 7. In the 
case of the non-precessing magnetization, the 7-modc 
corresponds to the longitudinal NMR mode. 



3. PARAMETRIC INSTABILITY OF HPD TO 
RADIATION OF 7-MODE 

31 Lagrangian for 7 mode 

In the simplest description, the dynamics of the 7- 
modc is determined by the following Lagrangian: 



C 



1 



X ( 7 2 -c 2 (V 7 ) 2 )+^(7) 



(12) 



Here we used the approximation of a single speed of spin 
waves c, since the effect of the anisotropy of the spin 
wave velocity has already been discussed in Ref. pQ . In 
the time-dependent part of Fd we only consider the first 
harmonic, i.e. according to Eq.Q we discuss the para- 
metric excitation of two 7- modes with ck w wl/2. The 
amplitude of the first harmonic is: 



Fx = — x^lsin/?! sin/3 2 cos( 7 /2) x 
2sl - 1 + ^ - ~ ~ + (1 + s)(l + I) COS7 ) . (13) 



Further we assume that the system is in the minimum of 
the dipole energy Fq as a function of 7. The equilibrium 
value 7 = 70 is 



cos 70 



(2aZ-l) + 2(l-a)(l-i) 
(l + s)(l + l) 



(14) 



which is valid if the right hand side of Eq. I|14|) does not 
exceed unity, i.e. when s + I — 5sl < 2. 

For the discussion of Suhl instability we need the 
time-dependent term which is quadratic in 7 — 70. Then 
the Lagrangian l|12|l which describes the parametric in- 
stability towards decay of Larmor precession to two 7- 
modes with kc ss wl/2 is (after the shift 7 — 70 — ► 7; 
neglecting Ql compared to u>l; and neglecting the 
anisotropy of the spin- wave velocity): 

C= i X (-7 2 + c 2 (V 7 ) 2 + aft 2 7 2 cos^) , (15) 
where, if s + I — 5sl < 2, the parameter a is 



a = — sin Pi sin fa 
15 



3(s + 1- si) 



2(l + «)(l + 
35, 



1/2 



(1 + s)(l + /) + 2(2sl - 1) + — (1 - s)(l - I) 



(16) 



32 Parametric instability 

Let us rewrite the Lagrangian (|15fl in terms of Hamil- 
tonian as function of creation and annihilation operators 
&k and 6£: 



Ok 



1 (bk ~ S) , Ws 2 (fc) = c 2 fc 2 , 



yj2xu s (k) 
Pk = X7k = v / X^(fc)/2(6 k + K) , 



H = J2 w.(*)&k&k + ( e_l " Lt& kfo-k + c.c.) 

k k 

V(k) = V Wf {k) = o^ L /u a (k) . 
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Here we neglected fi^ compared to u>l. The spectrum 
of the excited mode is b^(i) = b\ ti exp(—iu)i J /2t + v^t), 
where the instability increment v-^ is 

= VV(k) 2 -[io s (k)-uj L /2] 2 . (21) 

Increment reaches its maximum max v% = V{k) for 
spin waves, satisfying condition of the parametric 
resonance, uj s {k') = lu^/2. These modes grow exponen- 
tially: 

6 k (t) cx exp(V(k')t) . (22) 

At finite temperatures this growing is damped by 
dissipation, but at low temperature the dissipation be- 
comes small and catastrophic relaxation occurs. We 
shall follow Ref. and assume the spin diffusion mech- 
anism of dissipation. In this case the equation for tem- 
perature below which the instability of the homo- 
geneous precession towards radiation of spin waves with 
oj s (k') — ck' — ujl/2 starts to develop is: 

D(T cat ) = 2V(k')c 2 /Lu 2 L . (23) 

Here D(T) is the spin diffusion coefficient, which de- 
pends on temperature and decreases with decreasing T. 
In our case 

V BLV (fe') = 2an|/ Wi . (24) 

4. TWO MECHANISM OF SUHL INSTABILITY 

The characteristic feature of spin-orbit interaction in 
the isotropic 3 He-B is that it is symmetric with respect 
to the interchange of spin and orbital momenta 18 . 
This is the reason why Eqs. JHJ, (|13|) and Ijltjfl are 
symmetric under transformation (3\ <-> /?2 (or s I). 
The immediate consequence of this symmetry is that 
the Brinkman-Smith (BS) precession (/3i = (3, fa = 0) 
is equivalent to the static state of 3 He-B in magnetic 
field (fa = 0, fa = 0)- Since the static state is stable, 
the BS precession cannot be destabilized by spin-orbit 
interaction, which is also seen from Ea. (|16fl : the inter- 
action amplitude V(k) = if fa = (I = 1). 

This Z2 symmetry of the Larmor precession in 3 Hc- 
B is violated by the term in the gradient energy which is 
responsible for the dependence of the spin-wave velocity 
on the direction of propagation. This term, which we 
omitted in our consideration, becomes time dependent 
even in the background of the BS mode. This leads to 
the Sourovtsev-Fomin (SF) contribution V SF (k) to the 
interaction amplitude V(k), discussed by £Q for L||H. 
SF mechanism can be easily extended to the general 
precession with arbitrary fa. For that one must take 



into account, that the axis of anisotropy of spin- wave 
propagation is determined by the orbital vector L. This 
gives the following modification of V SF (k'): 

t/ ru>-\-^v tu'\ A^l sin/3i|l-2s| 

V 3F (k ) => V SF (k ) = ^— 2g2 _ 2g + 5 I sin 2S\ , (25) 

cos J = k-l . (26) 

Here \i — 1 — Cj_/c| is the anisotropy of the spin- wave ve- 
locity, where c\\ (c±) is the velocity of spin waves propa- 
gating along (perpendicular to) 1. We assume that both 
£1 2 l /lo 2 l and the spin- wave anisotropy /1 are small pa- 
rameters. In general 

V(k) = V BW (k) + V SF (k) . (27) 

Comparing these two contributions to V(k), causing 
the catastrophic relaxation - due to spin-orbit interac- 
tion, V(k) — V BLV (k), Eq. (|2*0|) and due to spin- wave 
anisotropy in Eq. (|2*5|l - one can see that V SF (k') > 
V BLV (fc') at high fields, when the Larmor frequency ujl 
essentially exceeds the Leggett frequency ^l- However, 
in the typical experiments |1U1 lll| the spin-orbit am- 
plitude, Vg LV (fc), appears to be more important. There 
are two reasons for that: the field is not sufficiently 
high (Ql/lol ~ 0.3 — 0.5); the texture of the orbital 1 
is typically formed due to boundary conditions at the 
walls of container: 1 must be oriented along the nor- 
mal to the wall. The ID simulation JJ] of the texture 
of the s and 1 fields (angles /3i and P2) under condi- 
tions of experiments but in the parallel-plate geometry 
is shown in Fig.l. The main part of the deviation from 
the BS mode develops in about 1 mm size region near 
the walls, which means that in the cylindrical cell the 
texture should be very similar. Thus the real mode of 
precession in the finite vessel is very different from the 
ideal BS configuration, which is static in the precessing 
frame. On the contrary, the angle 9 of the order param- 
eter R a i is oscillating as shown in the upper part of the 
Figure. These oscillations lead to parametric instability 
at low temperature. 

Figure 1 also shows the parameter a(x) as function 
of the coordinate x. Since the diffusion damping of spin 
waves occurs in the whole volume of the cell, the incre- 
ment V(fc') in Eq. (|24fl must be averaged over the cylin- 
drical cell: V(fc') = 2aVt 2 L /uo L , with a= 0.099 for our 
texture. This leads to the critical value of spin diffusion 
at which the catastrophic relaxation must occur 

ft 2 c 2 

D(T cat ) = 4a— \i = 0.085cm 2 /s (28) 

For £Il =244 kHz and ujl — 460kHz this value corre- 
sponds to the temperature of the catastrophic relaxation 
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Teat = 0.5T C , according to ^!5]i m a g°°d agreement 
with the experimental value T cat = 0.47T C ([HB IH])- 

In Ref. the critical diffusion for optimal con- 
figuration for spin wave velocity anisotropy was esti- 
mated as about 0.03 cm 2 /c. Because the texture of 
the 1- vector influences V SF (k') in Ea. (|25|) . the increment 
must be even smaller. This indicates that the spin- 
orbit mechanism of catastrophic relaxation is dominat- 
ing under the conditions of the experiments (101 111) . 
V BLV (k')>V SF (k>). 

To verify that it is really so, let us analyze the Greno- 
ble experiment 15 , in which the texture was destroyed 
by the RF pulse so that L was parallel to H even in 
the surface layer. Under these conditions the catas- 
trophic relaxation was observed at lower temperature 
of about 0.4 T c . In this configuration one has (3% = 0, 
the spin-orbit term, U BLV (fc), is switched off, the spin- 
wave anisotropy of Ref. [I], leading to V SF (k), becomes 
the main source of instability, and indeed this T cat well 
corresponds to the critical spin diffusion D(T ca ±) about 
0.03 cm 2 /c, calculated in Ref. [Jj. When the texture 
is restored, then if the spin-wave anisotropy is the only 
mechanism of the catastrophic relaxation, the texture 
should lead to decrease of V(k) and thus to decreasing 
T ca t. Instead, T cat increases demonstrating that when 
the texture appears another mechanism of instability 
emerges as we discussed here. 

Two contributions, Eq. (J37J, to the increment of 
parametric instability (|21[1 can be also compared using 
results of experiments made in Cornell In these 

experiments T cat was measured for 31 bar at different 



u)l. They found that T cat decreases from 0.39 T c at 
lo l = 600 kHz to 0.24 T c at uj l = 3 MHz. This is well 
described by the spin-orbit mechanism of Suhl instabil- 
ity, which becomes weaker at higher field according to 
Eq. i|22|) . On the contrary, if the spin- wave anisotropy 
is the dominating mechanism, T cat would not depend 
on u>l (there are processes in the spin- wave anisotropy 
mechanism whose increment is proportional to V1 2 l /uj\ 
Ref.pP as in the case of the spin-orbit mechanism, but 
they are relatively small, since contain the product of 
two small parameters, \i and Q? L /wfj . 



5. CONCLUSION 

We calculated in a unified way two contributions 
7 BLV (k) and V SF (fc) to the increment of the parametric 
instability of the Larmor precession in superfluid 3 He-B 
for any angle between L and H. The term V SF (k), pre- 
viously discussed in Ref. pp for simple geometry with 
L| |H, originates from the anisotropy of spin- wave veloc- 
ity. We modified this contribution, V SF (k) => Vg P (k) for 
the case of the texture of the order parameter, which 
occurs in real experiments. In addition we found an- 
other term, V BUV (k), which originates from the spin- 
orbit interaction. This second mechanism only takes 
place when the orbital momentum L deviates from its 
symmetric orientation along the magnetic field, in par- 
ticular in the presence of texture. Our analytical result 
for the onset of the parametric instability due to this 
mechanism is in a good quantitative agreement with 
experimental results. In particular it gives the correct 
dependence of the catastrophic relaxation temperature 
on magnetic field. The spin- wave anisotropy mechanism 
must dominate either when the texture is destroyed or 
in the limit of strong magnetic field. 
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Recherche", France. The work was also supported in 
part by the Russian Foundation for Fundamental Re- 
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